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Abstract 

An investigation of the processes Z — > gg and Z —* 77 in a background magnetic field is 

o 

presented. For homogeneous fields corrections to the charged fermion propagator can be calculated 
kjfy in leading orders of the magnetic field. This work examines the first order contributions of the 
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corrected propagator to decays that are otherwise zero. Results of the decay rates for varying field 



strengths are included. 
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I. INTRODUCTION 



The interaction of matter with magnetic fields has long been studied for the important 
consequences it has for particle dynamics. Whether it is from the astrophysical implication 
of a magnetic field on the structure of a neutron star or to the Zeeman shift of energy levels 
within an atom, the presence of an external magnetic field can open the door to interesting 
effects not previously seen. 

The effect of an homogeneous background field on the charged fermion propagator in QED 
was first studied by Schwinger Pjj. Schwinger's proper-time method provides the machinery 
for an exact solution to the altered propagator. Recently, a weak field expansion of this 
method has been applied to demonstrate that there exists enhancements to neutrino-photon 
interactions that involve charged fermion loops || [|, || . The goal of this article is to present 
what effects a magnetic field has on the decays of the Z which are otherwise forbidden. 

The Landau- Yang theorem uses Bose symmetry and rotational invariance arguments to 
prove that a vector particle cannot decay into two massless vector particles |5|, |J. That is, 
the decays Z — > 77 and Z — > gg are not allowed. However we can employ the field expansion 
of Schwinger's method to calculate if and to what degree the magnetic field stimulates the 
decay. 

We expand just to first order in the magnetic field B. This turns out to be a very good 
approximation considering that the relevant expansion parameter is \Qe\B/m 2 , where m is 
the mass of the fermion and Qe is its charge. For the lightest of the possible fermions, the 
electron, the critical field turns out to be of the order B = m? e /e ~ 4.4 • 10 9 T. This field 
is a great deal larger than those magnetic fields produced by very dense white dwarfs, and 
near the neutron star magnetic field intensity. The calculation is complete to first order, but 
because of this expansion parameter it is shown that the dominant contributions are from 
the lightest quarks (u, d) for Z —> gg and from the electron in the Z — ► 77 case. 

First we briefly introduce the reader to the methodology used in determining the changes 
in the fermion propagator. Next, the result of that technique is presented. This result is 
then applied to the decay mode of the Z involving two gluons. The results of this calculation 
are presented and contrasted with results from the similar decay into two photons. 
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II. OVERVIEW OF SCHWINGER'S RESULT 



Schwinger's method uses an integration over a proper-time variable to solve for the 
Green's function of the inhomogenous differential equation for a Dirac field in the pres- 
ence of the electromagnetic gauge field 

(i@-eQA'-m)G{x,x') = 5\x-x'). (1) 

G(x,x') is considered the matrix element of an operator G, such that G(x,x') = (x\G\x'). 
Therefore Eq. (D) is written as an operator equation with its indices suppressed 

(H-m)G=l, (2) 

or 

1 yi + m 



G 



fL — m H 2 — m 2 
—i 



/ ds(H + m)exp[i(H 2 - m 2 )s}. (3) 
Jo 

s is the proper-time variable, and 11^ is a generalized momentum operator defined by 

n„ = v» - eQA„. (4) 

This generalized momentum operator is characterized by the following commutation prop- 
erties 

[U^x v \ = ig^ (5a) 
[n^IL,] = -ieQF^, (5b) 

where is the electromagnetic field-strength tensor. 

The quantity exp[iH 2 s] is regarded as an operator that evolves the system in accordance 
with the "Hamiltonian" 7i = — J] 2 . That is, the matrix element of the operator U(s) = 
exp[i ~[l 2 s] is the time development function that takes the state x^(0) = x'^ to the state 

(x\U(s)\x'} = (x(s)\x(0)}. (6) 

This view of the transformation as an evolution of the state governed by a Hamiltonian 
parallels the idea that the fermion moves from one space-time location to another in a way 



that is determined by this Hamiltonian. Accordingly, there exist "equations of motion" 
satisfied by the generalized position and momentum operators 

dx a 



ds 
ds 



-i[U^,H] = 2QeF M jr. 



(7a) 
(7b) 



The relationship in Eq. (|7b|) only holds when F^ v is a constant. This is precisely the case 
in which we are concerned. For calculation purposes the magnetic field is considered to be 
constant and directed along the z-axis (F 12 = —F 2 \ = B, and F^ — for i,j ^ 1, 2). In 
addition there are three differential equations that describe the transformation function 



id a (x(s)\x(0)) = (x(s)\H\x(0)) 
{id^ + eQA^(x(s)\x(0)) = (a;(s)|n M (s)la;(0)> 
-^ + eQA M )(x( S )|x(0)> = (x( S )|n M (0)|x(0)>, 



(8a) 
(8b) 
(8c) 



and a boundary condition 



{x(s)\x(0))\^ = 5\x-x'). 



(9) 



Using Eqs. (|7]), (Bp, and (0) one can solve for the explicit form of the Green's function. 
For the purposes of this article, only the result will be given [|l(]]. Choosing a straight line 
to connect the space-time points x and x', the Green's function is 



G(x,x') = $(x,x') 



d 4 p 



-ip(x—x') 



Sip), 



(10) 



with 



$(x, x') = exp 



-leQ f 'd£M M (fl 

Jx' 



(11) 



and 



Qip) 



-i I dsexp 
'o 



2 2 tan(eQBs) 2 
—ts\ m — pi. H p ± 



eQBs 



(12) 



cos(eQBs) cos 2 (eQBs) 

The quantity p| indicates p% — p|, while pj_ = £>i + £>2- Similarly, ^ = 7 po — 73P3 and 
= JiPi + liVi- By expanding Eq. (|T2| ) to first order in the magnetic field the integration 



simplifies, and our result is 

By making use of the fact that the magnetic field is homogeneous, the result in Eq. ([13]) was 
cast into covariant notation. 



III. APPLICATION TO Z -» gg 

At this point the general results from the previous section are used to calculate the decay 
rate of the neutral Z boson into to two gluons. To lowest order in the perturbation theory 
there are two Feynman diagrams which contribute to this process; see Fig. |I| Let us first 
turn our attention to the leftmost contribution of Fig. |l[ 



* — -«T>fto75o^ k\ — <T>o75o7r> k\ 



k 2 ^*~^<S}QQQSL> k 2 



FIG. 1: The contributions to the decay of the Z with momentum k into two gluons with momenta 
k\ and hi- 

The T-Matrix for the first of our diagrams is written as 

r 3/2 



Ta = EE a a '^a tr(A"A b ) / d 4 x / d*y [ dS 
^^cos^sin^ V J J J y 

i a,b 

xtr [YG\x, y)YG\y, z)Y(g l v ~ aW)G\z, x)] 
xe^iK^K^exp [-i(kz - k lX - k 2 y)\. (14) 

The sum over i indicates the sum over all possible fermions that contribute to this loop 
diagram. The sums over a and b indicate the sums over the possible states of the gluons. 
However, one of these sums is eliminated by virtue of the trace over the two SU(3) generators 
tr(A a A 6 ) = 25 ab . From this point the summation over color will be represented by the 
repeated color index a. The couplings g l v and g\ are the vector and axial vector couplings, 
respectively; both are dependent on the particular fermion that couples to the Z boson. a s 
is the strong coupling constant evaluated at M z . 



Because there is the combination of the three Green's functions in Eq. fll4|) , there exists 



a product of three phases of the same form as in Eq. (|i~T|). As indicated previously, the 



integration from two space-time points is along a straight line. A convenient choice of gauge 
for this integration is one in which 



^4(0 = §(0,^2 -6, £1-^1,0). 



(15) 



Generally, the zeroth and third components of can be any constant but the above is 
chosen for its simplicity. With the choice of gauge found in Eq. fllS]) the total phase factor 
of the diagram is written 

ieQ 



y)&(y, z)&(z, x) = exp 



[z - yYF^z 



1 , ie Q ( 

1 +—z-{z 



(16) 



The fact that the magnetic field is constant and along the z-axis has been exploited. When 
the first order corrections to the Green's function found in Eqs. ( |T3D and fllBD are substituted 
into Eq. (|TJ]), the result is thought of as the sum of five separate contributions 

Ta = Tao + Ta P v1 + Ta P v2 + Ta P t3 + TAph- (17) 

There is the zeroth order contribution to the diagram Tao which is completely independent 
of the magnetic field. There is a first order correction to each of the legs of the fermion loop 
Ta P ti, Ta P t2 and Ta PT 3- And lastly, there is a first order correction to the overall phase of the 
diagram Ta p \i- Because the invariant amplitude of a given diagram M. is defined through 
the relationship T = {2-r[) A AA5 A (k — k\ — k 2 ), it is straightforward to show that the zeroth 
order contribution from the first diagram is 



M 



AQ 



cos 6w sin 6w 
d 4 p 



X 



(2tt)4 



tr - my l Y($ + fa - m)' 1 



m) 



(18) 



For the first order correction to one of the legs of the loop, the contribution to the amplitude 
is 

-i7i- 3 / 2 eF Q(3 a sy fa 



M 



A P rl 



cos 9w sin 6 



d 4 p 



tr 



r 



<rFi + m ) 

(p 2 — m 2 ) 2 



i v {i> + fa 



m) 



m) 



(19) 
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The corrections to the remaining legs of the fermion loop are similar to Eq. ([19]). However, 
for the contribution due to the phase it is necessary we recognize that 



T Ap h oc - yY(z - x) u e 



" -ip '{y-z) e -ip" '(z-x) 



xtr 



(20) 



where p' is the momentum carried by the fermion as it moves from point z to point y, and 
p" is its momentum from point x to z. The relationship in Eq. (|20|) is rewritten in terms of 
derivatives with respect to the fermion momenta 

d 8 



TAph. oc 



dp'^ dp" v 



e ~ip' (y-z) e -ip" (z-x) 



xtr 



r 



(21) 



By making use of integration by parts and the identity 

|c-. w = _ C - 1(I) ^) C - IW , 

the first order contribution of the correction to the phase to the amplitude is 

— i^^eFapasy/a ■ 



M 



Aph 



X 



cos 9w sin 6 W 
d 4 p 



(2vr) 4 
X7 ff (^-^7 5 



tr[7 M (^-m)~ 1 7 I '(/+ ^ 2 -m)~ 1 -f l3 (p / + ^ 2 - m) -1 



(22) 



We follow the same procedure for the second diagram of Fig. [I] as we did for the first. 
The results obtained are similar to those found in Eqs. (|18D, (|19"D, and (^) with only the 
labels (1, 2) interchanged. 



IV. THE DECAY RATE 



The sum of the zeroth order contributions from the two diagrams vanishes. This is 
consistent with the fact that the process Z — > gg is forbidden by the Landau- Yang Theorem 
in the absence of a magnetic field [[5], The use of the charge conjugation operator simplifies 
the contributions that are first order in the magnetic field. We find that the axial vector 
parts of the two diagrams cancel while the vector parts are identical and add. Tracing over 
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the Dirac indices and integrating over the loop momentum results in an amplitude that is 
proportional to forty-four separate tensor structures formed by possible contractions of the 
polarization vectors, the field-strength tensor, and the gluon momenta. Since the amplitude 
is symmetric under interchange of the labels (1, 2), we find that the tensor structures can be 
grouped into twenty-two pairs; each pair having its own unitless coefficient which depends 
on the possible fermion masses, charges, and couplings in the charge loop. Limiting the 
gluons to be only transversely polarized allows us to drop all terms proportional to (e* a • fcj) . 
This cuts the number of tensor structures in half to eleven. Thus, the amplitude is written 
as 

+N 2 [(e? ■ e)(e* 2 a Fh) + (e? ■ e)(el a Fk 2 )] (1/M Z ) 

+N 3 [(el a ■ e)(e?Fk 2 ) + (e? ■ e)(eTFh)] (1/M Z ) 

+N 4 (e* 1 a -e* 2 a )(eF(k 1 + k 2 ))(l/M z ) 

+N 5 (e? ■ e* a )(e ■ (h - k^ihFk^l/Mzf 

+N e (el a ■ k 2 )(e* a • fc 1 )(e • (h - k 2 ))(hFk 2 )(l/M z ) 5 

+N 7 [(e{ a ■ e)(e* 2 a ■ h) - (e? ■ e)(e? ■ k 2 )] (hFk 2 )(l/M z ) 3 

+N 8 [(e* a ■ h)(e ■ k 2 )(e{ a Fk 2 ) + (ef • k 2 )(e ■ k^Fk^l/Mzf 

+N 9 [(e* 2 a ■ h)(e ■ h)(el a Fk 2 ) + (ef • k 2 )(e ■ k^Fk^l/Mzf 

+N 10 [(e* a ■ h^e ■ k^Fh) + (ej a • k 2 )(e ■ h)(e* 2 a Fk 2 )] (1/M z f 

+N u (el a ■ k 2 )(e* 2 a ■ fc 1 )(eF(fc 1 + k 2 )) (1/M Z ) 3 }, (23) 

where M z is the mass of the Z, 9 W is the Weinberg angle, and the iVj are the unitless 
coefficients. The reader should see the appendix for a detailed explanation of how the 
coefficients are defined. 

At this point we confirm that our amplitude is gauge invariant by replacing a polarization 
vector by its associated momentum 

eT{h) - K , (24) 

and verifying that the amplitude vanishes. For these purposes we drop the color index a. 
The gauge invariance implies six independent relations among the coefficients that reduce 
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the number of independent coefficients from eleven to five. These six relations and five 
coefficients are given explicitly in Eqs. ( |A1| ) and ( |A2] ) of the appendix. 

What is of real interest, however, is the square of this amplitude. When the square is 
taken and evaluated in the center of mass frame of the Z, we find that the square of the 
amplitude depends on only three of our five independent tensor coefficients 

\Mz^ 99 ? = — J a if ■ 2fl (4MWg+|iV a -iV a |Wg), (25) 
2567T MJ cos^ Ow sin 9 W 

where 9 is the polar angle as measured from the direction in which the magnetic field points. 

It is useful to extract from Eq. (|25|) the effective expansion parameter for the magnetic 
field. That is, it was discussed earlier that the perturbation theory for the magnetic field 
is an expansion in \Qe\B /m? . Because this is fermion dependent, it is useful to pick the 
lightest fermion and make B/Bq the expansion parameter. Recall that Bq = m? e /e. Making 
this substitution, Eq. (p5|) becomes 



\M^ s f = :f- A ff% (^)V,|W» + |AT, - ^We). (26) 
2567T cos z 9w sin Ow \M Z J 

At this point it is instructive to note that the leading contributions of N\ go as 
while for both N2 and N% the leading contributions go as 

(^f) E3- Qualitatively, this 
means that only the lightest of the fermions make significant contributions to the the decay 
rate. It also means that the Z is many times more likely to decay such that the gluons 
travel perpendicular to the magnetic field rather than parallel. Computation reveals that 
this factor is 

2 . io 13 > \ N2 ~ N3 \ 2 > q . 10 io 
U ~ 4|JVi|2 ~ D iU ' 

The large range in values is a symptom of the very large uncertainty in the current quark 
masses. If the constituent quark masses are used, the ratio drops to a factor of 2 • 10 4 . 
The decay rate for this process is 

r z _ = r^y'' te)W + liW,l'). (27) 

3 • 2 n 7i 2 cos 2 9 w sm 2 9w \M Z J K ' 

Because of the large current quark mass ranges quoted by the Particle Data Group 0], there 
is a corresponding large range in possible values for the decay rate. We choose to use the 
mean values of these masses for computation. Seperate calculations are carried out using 
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the constituent quark masses. Given below are the results of the decay rate for both the 
current and constituent mass calculations fll2"| , respectively, 



9.60 ■ 10" 15 \ B < 




. 5.10 -10- 7 . , ., . 

r ^,= _ 1E hr MeV - ( 28 ) 



0, 



The values of the branching ratio of the above decay rate with the total width of the Z in 
the absence of a magnetic field are 



Z->gg 



2.04 • 10~ 10 
3.85 ■ 10~ 18 




(29) 



For comparison purposes it is convenient to contrast the results of the two gluon case 
with the two photon case. The results of the latter decay were determined in a manner 
completely analogous to that which has been outlined in this article for the decay Z — > gg. 
The values of the decay rate for the process Z — > 77 are 

. 8.33 -lO^ 8 . , „ . 
Y z ^ n ={ o W — MeV. (30) 



6.63-10 



-8 




Notice that the degree to which the two values differ is much smaller than in the two gluon 
case. This is because the major contribution to the process is due to the fermion mass being 
that of the electron, which is known very well. The corresponding values for the branching 
ratio are 



Z-»77 



3.34- 10~ n 
2.66 • 10~ n 




(31) 



V. CONCLUSION 

This article has demonstrated that the existence of a magnetic field can stimulate de- 
cays that were otherwise forbidden. Schwinger's powerful proper-time method was used in 
conjunction with a first order expansion in the field to determine what corrections exist for 
the Dirac Green's function. These new contributions were analyzed in the specific decays 
Z — > gg and Z — > 77. Though the resulting decay rates are small, they are non-zero. In 
addition it was also found that there exists a very large asymmetry as to whether the decay 
products will exit perpendicular or parallel to the magnetic field. 
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APPENDIX A 

As noted in the body of the article, the amplitude for the decay Z — > gg to first order 
in an external magnetic field (Fig. |l]) is the sum the contributions from the zeroth order 
amplitude (Eq. fll8|)), first order corrections to each leg of the charged ferion loop (Eq. ([19]) 
and the like), the first order correction to the phase of the amplitude (Eq. (p2|)), and similar 
corrections for the crossed process which one obtains by interchange of the labels (1,2). 

Evaluation of this amplitude gives the result found in Eq. (P3|). Gauge invariance of this 
amplitude implies the following six relationships which we confirmed, 



Nq = AN Y - 4A 2 - 2N 5 (Ala) 

N 7 = 2N 2 (Alb) 

N 8 = 2N 1 -2N 2 (Ale) 

N 9 = -2N l (Aid) 

N 10 = -2N 3 (Ale) 

A n = -2A 4 . (Alf) 
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Explicitly, the five remaining independent coefficients are 



r ( ( < - 4 )' 



Qigi 



Qi9l 

(*-4) J 



(*-4)' 



(-128*? + 768^- 1024) 

+ (96t- - 896 U + 2048) 
- (64 ^ - 512 + 1024^ 1 )S(^) 

(-8 tl + 32 t\ - 256 ti + 1024) 

-(80^-448^ + 512)^) 

-(l6t 4 2 - 64 ti - 256 + 1024 t^Bfa) 

(-(8/3) if + (160/3) t f ?- (512/3) *0 

+ (48t 2 - 320t; + 512) 
-(16^-128^ + 256)5^) 

(8^-64^ + 128^) 

-(48<f-384*i + 768)A(ti) 
-(16^-128^ + 256)5^) 

(16 t\ + 64 1 2 - 1024 ^ + 2048) 

-(224 - 1920 U + 4096) A(^) 

+ (32t 2 - 128 ti - 512 + 204&ti 1 )B(t i ) 



where ti is defined by the square of the ratio of the Z mass to the quark mass 



Mo 



mi 



and the functions A(ti) and B(ti) are defined by 



A(t) = j^dx ln[l - tx{l -x)-i8] = 2 



'^sinh" 1 J=£-\ 



^ sin" 1 J I - 1 



(A2a) 



(A2b) 



(A2c) 



(A2d) 



(A2e) 



(A3) 



cosh 



-1 It ITT 

4 2 



£ - 1 



for t < 



for 4 > t > 



, for t > 4, 



(A4) 



12 



and 



sinh 1 \l 



sin 



-1 / i 



for t < 



for 4 > i > (A5) 



cosh 



1 2 



1 t 



7T 2 - Am cosh \/|, for t > 4. 



Notice that in Eqs. ( |A2|) there is a sum over all possible quarks in the charge loop. For 
the decay Z — > 77 this could also include the electron. and g l v are the charge and the 
associated vector coupling of the fermion, repectively. The factors of Mz present are due to 
the evaluation of the dot product between the two gluon momenta 



ki ■ k 2 



(h + k 2 f/2 

k 2 /2 

Mf/2, 



(A6) 



where k is the four-momentum of the Z. 
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